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ABSTRACT 



Pathological folding and oligomer formation of the amyloid /3-protein (A/3) are widely per- 
ceived as central to Alzheimer's disease (AD). Experimental approaches to study A/3 self- 
assembly provide limited information because most relevant aggregates are quasi-stable and 
inhomogeneous. We apply a discrete molecular dynamics (DMD) approach combined with 
a four-bead protein model to study oligomer formation of A/3. We address the differences 
between the two most common A/3 alloforms, A/340 and A/342, which oligomerize differ- 
ently in vitro. Our previous study showed that, despite simplifications, our DMD approach 
accounts for the experimentally observed differences between A/340 and A/342 and yields 
structural predictions amenable to in vitro testing. Here we study how the presence of elec- 
trostatic interactions (Els) between pairs of charged amino acids affects A/340 and A/342 
oligomer formation. Our results indicate that Els promote formation of larger oligomers in 
both A/340 and A/342. Both A/340 and A/342 display a peak at trimers/tetramers, but A/342 
displays additional peaks at nonamers and tetradecamers. Els thus shift the oligomer size 
distributions to larger oligomers. Nonetheless, the A/340 size distribution remains unimodal, 
whereas the A/342 distribution is trimodal, as observed experimentally. We show that struc- 
tural differences between A/340 and A/342 that already appear in the monomer folding, are 
not affected by Els. A/342 folded structure is characterized by a turn in the C-terminus that 
is not present in A/340. We show that the same C-terminal region is also responsible for 
the strongest intermolecular contacts in A/342 pentamers and larger oligomers. Our results 
suggest that this C-terminal region plays a key role in the formation of A/342 oligomers and 
the relative importance of this region increases in the presence of Els. These results suggest 
that inhibitors targeting the C-terminal region of A/342 oligomers may be able to prevent 
oligomer formation or structurally modify the assemblies to reduce their toxicity. 
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Introduction 



Alzheimer's disease (AD) is a progressive brain disorder, clinically characterized by the 
accumulation of extracellular amyloid deposits composed of amyloid /5-protein (A/3), intra- 
cellular neurofibrillary tangles, and neuronal loss. Recent research supports the hypothesis 
that cerebral A/5 accumulation is the primary cause of neurotoxicity in AD Accumulating 
evidence suggests that A/3 oligomers and prefibrillar aggregates are the proximal effectors of 
neurotoxicity in the early stages of AD . The predominant forms of A/5 found in brains 
of AD patients are 40 amino acids long (A/340) and 42 amino acids long (A/542). A/542 is 
linked particularly strongly with AD. Genetic studies have shown that autosomal dominant 
forms of AD invariably involve increased production of A/5 or an increased A/542/ A/340 con- 
centration ratio |4|. A/542 forms fibrils at significantly higher rates than does A/540 

BB 

and A/342 self-association produces structures that are more neurotoxic than those formed 
by A/540 3j. Experimentally, there is a distinct difference in oligomerization pathways of 
A/540 and A/542 2j- In vitro experiment using the techniques, photo-induced cross-linking of 
unmodified proteins (PICUP), size-exclusion chromatography, dynamic light scattering, cir- 
cular dichroism spectroscopy, and electron microscopy showed that A/5 exists as monomers, 
dimers, trimers, tetramers, and larger oligomers in rapid equilibrium. The A/540 oligomer 
size distribution comprises monomer, dimer, trimer, and tetramer, in similar amounts, and 
few higher-order oligomers. The A/542 distribution is multimodal, displaying a prominent 
peak of pentamers/hexamers and smaller peaks of dodecamers and octadecamers 0]. 

Detailed, quantitative analysis of the three-dimensional structures, energetics, and dy- 
namics of oligomer formation are necessary steps toward a molecular understanding of A/5 
assembly and neurotoxicity. During the formation of fibrils, oligomers of different sizes co- 
exist with monomers and larger aggregates such as protofibrils [8J and fibrils. The relative 
amounts of each oligomer type are small, which makes determination of the structural prop- 
erties of the oligomers difficult. Computer simulations, in contrast, are not subject to the 
same kinds of problems, allowing small oligomers to be studied at full atomic resolution (for 
recent reviews, see 10 1, and 11]). 

Conventional "all-atom" molecular dynamics (all-atom MD) simulations with explicit 
solvent which take account of all the protein and solvent atoms give the most detailed 
information. However, aggregation studies using all-atom MD with explicit solvent are 
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currently limited to either aggregation of small number of A/3 fragments such as three A/3 (16 



13 



22) peptides [12j or stability studies of various A/3 dimers with predetermined structures 
\h\ . Tarus et al. used a protocol based on shape complementarity to determine the initial 
A/3io_35 dimer structure and showed that the peptide dimers are stabilized primarily through 
hydrophobic interactions 13j. Huet et al. studied A/340 and A/342 dimers and their A21G 
conformers, starting from their fibrillar conformations and found various possible topologies 
of dimers in equilibrium [l^] . Keeping track of positions and velocities of all the atoms at 



the times simulated by the 



101 ] . However, protein folding 



every time step is computationally expensive. Consequent 
all-atom MD simulations are limited to few microseconds 
and aggregation usually occur on time scales larger than milliseconds. To overcome this 
limitation, we use fast and efficient discrete molecular dynamics (DMD) simulations 3] 
with a simplified four-bead protein model and implicit solvent. DMD is a simplified version 
of MD using combination of square-well potentials. The DMD approach with a simplified 



protein folding and aggregation 



protein model and implicit solvent increases the efficiency o 
studies by a factor of ~10 7 compared to the all-atom MD 161 ] . 

The idea of app lying the DMD approach to study protein folding was proposed in 1996 



by Zhou et al. 



171 ]. Soon after, the method was combined with a one-bead protein model to 



18 



19 



20 



21 



22j. In 2004, Peng et al. 



study folding of a model three-helix bundle protein 
used DMD with two-bead protein model to study aggregation of an ensemble of 28 A/340 
peptides into a fibrillar structure 23|]. Smith and Hall introduced four-bead protein model 
in combination with the DMD, and showed a cooperative transition of a polyalanine chain 
into an a-helical conformation without any a priori knowledge of the native state 24, 125 ]. 
Using the four-bead protein model and hydrogen bond interactions in combination with the 
DMD on a single 16-residue polyalanine chain, Ding et al. demonstrated a temperature- 
induced conformational change from the a-helix to the /3-hairpin conformation 26J. Urbane 
et al. studied folding and dimer formation using DMD with the four-bead protein model, 
and investigated stability of dimer conformations predicted by DMD approach using all- 



atom MD simulations 



271 ] . Lam et al. used the same model to study the A/342 folding 



and its temperature dependence 



2s 



The results of Lam et al. were in a good qualitative 
agreement with an all-atom study using implicit solvent 29j and, importantly, consistent 
with the temperature dependence of A secondary structure, experimentally determined by 
Gursky and Aleshkov and Lim et al. 30, \3l\\ . 
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Recently, we studied oligomer formation using a four-bead model with backbone hydro- 



gen 
Els 



jond interactions and the amino acid-specific hydropathic interactions, but no effective 



321 ] . We observed that dimers are the most abundant among the low molecular weight 



A/540 oligomers and that the frequency of trimers and higher-order oligomers decreases 
monotonically. In contrast, the A/542 oligomer size distribution was bimodal, with signifi- 
cantly more pentamers than A/540. Multimodal and unimodal oligomer size distributions are 
discriminating properties of A/542 and A/540, respectively, as observed in vitro by PICUP Q. 
Experimentally-detected pentamer/hexamer A/542 oligomers were termed paranuclei. Ex- 
istence of A/542 paranuclei and their homotypical assemblies, "oligo-paranuclei" , has been 
independently confirmed by a combination of ion mobility and mass spectrometry {33]. Im- 
portantly, paranucleus-like assemblies have been detected in vivo in the form of dodecameric 
assemblies termed ADDLs 



34] . globulomers |35|, and A/5*56 



that oxidation of M35 blocks A342 paranucleus formation 37] . A/3 without oxidated M35 



In vitro studies showed 



displays both accelerated [38|, |39[ and delayed [40] fibrillogenesis rate relative to wild type 



A/5. Analysis of intramolecular contacts in A/540 and A/542 pentamers in our in silico study 



also showed that M35 forms contacts with 141 and A42 in A/542 
nation of the above experimental results 



321 ] . providing an expla- 



371 ] . In addition, our prior study indicated that 



A/542 monomers but not A/540 monomers are characterized by a turn structure, centered 
at G37-G38, and that this turn structure was more prominent in large oligomers 32]. This 
result is consistent with recent proteolysis results using A/540 and A/542 41]. 



There is indirect in vitro as well as in silico evidence suggestin g th at 



role in both A/5 folding 41 



42 



43 



44] and A/5 fibril formation 45|, 



46 



Is play a significant 



47] . In the present 



study, we follow the protocols of our previous study [32J using DMD and four-bead protein 



model with amino acid-specific interactions 



ll] to elucidate the role of Els between pairs of 



charged amino acids (D, E, K, and R) on folding and oligomerization of A/540 and A/542. 



Methods 

For our simulation method, we use DMD simulations [lj]] . The main simplification in this 
method is to replace continuous interparticle potentials by a square-well or a combination 
of square- well potentials. As a result, particles move along straight lines with constant 
velocities until a pair of particles reaches a distance at which the interparticle potential is 
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discontinuous. A collision event then takes place during which the velocities and directions 
of the particles are updated while preserving the total kinetic energy, momenta, and angular 
momenta. Because DMD is event-driven, it is faster than all-atom MD. Our DMD approach 



111- 



using coarse-grained protein models has been described in detail elsewhere 

Here we use a four-bead protein model with backbone hydrogen bonding, effective hy- 
dropathic interactions and Els. We use the four-bead model with hydrogen bonding, intro- 



duced by Ding et al. 



261 ]. then further generalized by Urbane et al. 32| to include amino 



acid-specific hydropathic and electrostatic interactions. In the four-bead model, the back- 
bone is represented by three beads, corresponding to the amide (N), the a-carbon (C a ), 
and the carbonyl (C) groups. Each side-chain is represented by one bead (Cp). G, which 
lacks a side-chain, has no Cp bead. As the carbonyl oxygen and the amide hydrogen are not 
explicitly present, an effective backbone hydrogen bond is introduced between the nitrogen 
atom Ni of the i — th amino acid and the carbon atom Cj of the j — th amino acid. Because 
the solvent is not explicitly present in our DMD approach, effective interactions between 
the side-chain atoms are introduced to mimic the solvent effects. The relative strength of 
hydropathic interactions between pairs of side chain beads is based on the Kyte-Doolittle 



hydropathy scale 



48|. When two hydrophobic side chain beads are within the interaction 



range of 0.75 nm, they interact through a one-step attractive potential. When two hy- 
drophilic side chain beads are within the same interaction distance, they interact through a 
one-step repulsive potential. In our model, the hydrophobic amino acids are A, C, F, L, M, 
I, and V. The hydrophilic amino acids are D, E, H, K, N, Q, and R. The side chains of the 
remaining amino acids G, P, S, T, W, and Y interact only through a hard-core repulsion. 
The Els are implemented by assigning a two-step potential with two interaction distances, 



0.60 nm and 0.75 nm, as described elsewhere 



ll| . When two beads with the same charge are 



at the interaction distance, they interact through a positive (repulsive) two-step potential. 
Two oppositely charged beads interact through a negative (attractive) two-step potential. 

We set the potential e nerg y of the hydrogen bond, Ehb, which in proteins is typically in 
the range 1 — 5 kcal/mol 49(, to unit energy (Ehb — !)■ We set the potential energy of the 
hydropathic interactions Ehp = 0.3. Experimental free energy of salt bridge formation is 
estimated to be in the range 0.7 — 1.7 kcal/mol 50|, thus we choose the potential energy of 
Els, Eqh — 0.6. Using the unit of temperature EhbI^-b where &b is Boltzmann's constant, 
we estimate that T = 0.15 is appropriate for simulating physiological temperatures. We 
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perform DMD simulations in the canonical ensemble (NVT) using the Berendsen thermostat 



algorithm 



5l|- 



Because we treat the solvent in our DMD approach implicitly, the effective interactions 
between the side-chain beads include not only protein-protein but also protein-solvent in- 
teractions. Thus, there are no generic interaction parameters that would be independent 
of the environment. Moreover, because different proteins may interact with the solvent in 
different ways, the implicit effect of the solvent and thus the interaction parameters may 
depend on the particular protein sequence. The complexity of protein-protein and protein- 
solvent interactions represents a challenge in protein structure prediction where even the 
most successful specialized models fail on certain targets 52] . The question of how general 
is a particular choice of interaction parameters in our DMD approach is a topics of future 
studies. 



Results and Discussion 

We apply the four-bead model with hydrogen bonding and amino acid-specific interac- 
tions due to hydropathy and charge and use DMD with implicit solvent to study A/340 
and A/342 oligomer formation. Due to simplifications in protein description and implicit 
solvent, our DMD approach is efficient enough to allow for a study of the whole process 
starting from unfolded separated peptides to formation of quasi-stable A/3 oligomers with 
well-defined size distributions. In our protein model, each side chain is replaced by at most 
one bead, a significant simplification considering side-chain diversity. However, recent devel- 



opments in understanding of protein folding and assembly show that despi 



of the process as a whole, the underlying fundamental physics is simple 



53 



;e the complexity 



54| . It is believed 



that the patterns of hydrophobic and hydrophilic residues, rather than the highly specific 



characters of the individual residues involved, play an important role [55|, |56[]. This is consis- 
tent with our prior simulation results where we showed that amino acid-specific interactions 
due to hydropathy itself are sufficient 32] for accounting for the experimentally observed 



oligomer size distribution differences between A/340 and A/542. Here, we apply the same 
model, with the addition of Coulombic interactions between pairs of charged amino acids, 
to study the effect of Els on A/340 and A/342 oligomer formation. 

The primary structure of A/342 is DAEFRHDSGYEVHHQKLVFFAEDVG 
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SNKGAIIGLMVGGVVIA. The primary structure of A/340 is identical, except that the last 
two amino acids, I and A, are missing. We define the following peptide regions: (i) the 
N-terminal region D1-K16 (NTR); (ii) the central hydrophobic cluster L17-A21 (CHC); (Hi) 
the turn A region E22-G29 (TRA); (iv) the mid-hydrophobic region A30-M35 (MHR); (v) 
the turn B region V36-V39 (TRB); and (in) the C-terminal region V40/V40-A42 (CTR). 
The CTR of A/340 consists of only one amino acid, V40. 

We simulate eight oligomerization trajectories for A/340 and A/342 each, starting from 
spatially separated peptides. Each initial configuration consists of 32 A/340 (A/342) peptides 
with a zero potential energy and with randomized spatial positions and randomized initial 
velocities of atoms within a cubic box of side 25 nm. The molar concentration is ~ 3.4 mM. 
This initial setup follows the protocol of our prior publication 



321 ] . The concentration in 
our simulation is 10 — 100 times higher than that studied experimentally Lowering the 
concentration is possible only at a high cost of efficiency of our approach. As shown in a 
recent study by Nguyen and Hall 57J, lowering the concentration may give rise to a-helical 
aggregates at low temperatures, possibly altering the assembly pathways, a problem to be 
addressed in future studies. 

The energy is in our approach normalized to the potential energy of the hydrogen 
bond Ehb = 1- Temperature is expressed in units of energy and also normalized to 
Ehb- The maximal potential energy of the hydrophobic/hydrophilic interaction is set to 
Ech = Q-Q/Ehb = 0.6. The N-terminal amine group and the C-terminal carboxyl group 
are non-charged. 

Hydrogen bonding is the same for all amino acids and represents the basic interaction 
needed to model the secondary structure, a-helix and /3-strand, formation. When only 
the hydrogen bond interactions are allowed (Ehb = ^,Ehp = 0, and Ech = 0), a single 



planar /3-sheet aggregate is formed 111, |27(]. Thus, only hydrogen bond interaction is not 



enough for description of spherical oligomers with only small amounts of secondary structure. 
Recently, we introduced the effective hydrophobic/hydrophilic interactions which are amino 
acid-specific to mimic the effect of aqueous solution [32j . Using the hydrogen bonding and 
effective hydropathic interactions but no Els (Ehb — 1, Ehp = 0.3, and Ech = 0), we found 
spherical A/3 aggregates with a dense hydrophobic core and with the hydrophilic N-termini 



comprising the surface 



32]. 



The aim of the present study is to explore the effects of Els on oligomer formation of 
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A/340 and A/342. The question of how Els affect the aggregation is intriguing because most 
of the charged amino acids are at the N-part of the molecule: six of nine charged amino 
acids are within the D1-K16 fragment as opposed to the hydrophobic residues which are 
concentrated in the remaining fragment L17-V40/A42. Fig. [T]shows typical conformations of 
a folded monomer, dimer, and pentamer of A/342 in the absence and presence of Els. Similar 
conformations are found in the case of A/340 (data not shown). We observe various topologies 



at a fixed oligomer size, which is consistent with findings by Huet et al. [14] . To gain more 
quantitative insight into the oligomer formation and structure, we quantify the oligomer 
size distributions, calculate the intra- and intermolecular contact maps, secondary structure 
propensities, and Ramachandran plots for each A/340 and A/342 alloform separately. 

A/340 and A/342 oligomer size distributions 

All simulations are 10 million simulation steps long. Initially, all the oligomer size distri- 
butions are peaked at monomers and the oligomer size distributions of A/340 and A/342 are 
equivalent. The difference between A/340 and A/342 size distributions develops steadily with 
increasing simulation time and at ~ 6 million steps the difference between A/340 and A/342 
oligomer size distributions becomes statistically significant as determined by applying the 



X 2 -test 58j. When comparing oligomer size distributions of each alloform separately at 8.0, 
8.5, 9.0, 9.5, and 10.0 million steps, we find that within this time window size distributions 
do not differ significantly. However, the number of monomers and oligomers of all sizes 
is variable. Each of the final oligomer size distributions is obtained by first average over 
all 8 trajectories at a fixed simulation time, and then the resulting ensemble averages are 
averaged over the simulation times of 8.0, 8.5, 9.0, 9.5, and 10.0 million steps. 

We have shown previously that A/340 and A/342 oligomer size distributions in the absence 



of Els (Ech = 0) are significantly different (Fig. 12(a)) [32fl. A/340 and A/342 oligomer size 
distributions in the presence of Els (Ech — 0.6) are significantly shifted towards larger 
oligomers, as shown in Fig. [2(b). Comparing the A/340 and A/342 oligomer size distributions 
by applying the x 2 -^est, we conclude that in the presence of Els, the distributions are 
significantly different (p < 0.01). 

In the presence of Els, the average size of A/340 oligomers increases from 3.0 to 5.2 
molecules, and the average size of A/342 oligomers increases from 3.7 to 6.2 molecules. These 
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results suggest that Els facilitate aggregation. A/342 forms significantly more nonamers and 
larger oligomers compared to A/340. The A/340 size distribution is unimodal with a peak at 
tetramers. The A/342 distribution contains a trimer peak and two additional peaks, at n = 9 
(nonamer) and n = 14 (tetradecamer) , neither of which is present in the A/340 distribution. 
A multimodal oligomer size distribution was observed experimentally with A/342, but not 
with A/340 J. 

In our simulations, the N- and C-termini are uncharged, whereas in the experimental 
studies, the N-terminus is positively charged (NH^) and the C-terminus is negatively 
charged (COO—) Q|. Observation of high-order oligomers in our simulations is con- 
sistent with in vitro results in which the C-terminal carboxyl group was replaced by the 
electrostatically neutral carboxamide, resulting in a greater abundance of high molecular 
weight oligomers j^]. Our simulation results, in combination with experimental findings, 
thus suggest that inclusion of charged termini, in particular the C-terminal negative charge, 
will moderate formation of A/3 oligomers. This hypothesis will be tested in future compu- 
tational and experimental studies. 



Secondary structure of A/3 monomers 



We calculate the secondary 



structure propensities on each folded monomer separately 



using the STRIDE program 59] and then average over different conformations to obtain the 
average values of the a-helix, turn, and /3-strand propensities per amino acid. At 1 million 
(M) step, the potential energy of individual monomers is stabilized (data not shown), thus 
we consider monomers to be in a folded state at 1M step. 

Folded monomers do not have a significant amount of a-helix structure (data not shown). 
Figs. [3] (a) and (b) show the turn propensity per amino acid for folded A/340 and A/342 
monomers in the absence and presence of Els. A dramatic effect of Els on the turn propen- 
sities in both alloforms is observed in the region A21-A30. In the absence of Els this region is 
characterized by two turns, the first at A21-V24 and the second at S26-G29. In the presence 
of Els, only a single turn within the region V24-G29 remains. 

Figs. [3] (c) and (d) show the /3-strand propensity per amino acid for folded A/340 and 
A/342 monomers in the absence and presence of Els. As a result of Els in both alloforms, the 
regions A21-D23 and K28-I31 show an increased /3-strand propensity. In A/340 monomers 
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the regions A2-F4 and L34-G38 show a decreased /3-strand propensity due to Els. In A/342 
monomers the regions R5-H6 and L34-V39 show a slightly decreased /3-strand propensity 
due to Els. Notice that the /3-strand propensity per amino acid is below 40% for A/340 and 
below 30% for A/542. The number of turns and consequently also the number of /3-strand 
regions in the A/342 monomer (5) is bigger than in the A/340 monomer (4), indicating a more 
compact structure of the A/342 monomer as compared to the A/340 monomer, a consequence 
of a strongly hydrophobic CTR in A/342, which introduces an additional turn centered at 
G37-G38. The average turn and /3-strand contents of A/340 and A/342 folded monomers 
are displayed in Table I. These contents are calculated from propensities per residue by 
averaging over all residues in the the peptide. Table I shows that for both A/340 and A/342 
the average turn content is in the range 43-45% while the average /3-strand content is in 
the range 10-12%. Neither the average turn nor the average /3-strand content is strongly 
affected by Els. 

The above results suggest that even in the presence of Els, the A/3 monomer is a collapsed 
coil with several turns and some /3-strand but no a-helical structure, which is in agreement 
with existing experimental studies 30, 41, 6j|. The /3-strand propensity of A/340 monomer 



as shown in Fig. 3(c) is also consistent with a recent study of A/340 folding using a scanning 



tunnelling microscopy that showed monomers folded into 3 or 4 domains with some /3-strand 



structure 



fill- 



Intramolecular contacts of folded A/3 monomers 

Here we discuss the effect of Els on the intramolecular contacts among pairs of amino acids 
of folded monomers. Initially, monomer peptides are in zero-potential energy (unfolded) 
conformations. At 0.1M steps, over 60% of peptides (65.9% for A/340, 60.5% for A/342) are 
folded. We describe the regions of the most important contacts between pairs of amino acids. 
We first describe "short-range" contacts formed within the regions TRB, MHR, and TRA. 
Then, we describe the "long-range" contacts between the regions CHC-CTR, CHC-MHR, 
and CHC-NTR. 



Previous results [32| showed that while A/340 and A/342 monomers both display strong 
contacts within the TRA region, strong contact in the TRB region with a turn centered at 
G37-G38 are characteristic of A/342 only. This in silico difference between A/340 and A/342 
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folding is consistent with experimental findings by Lazo et al. 41]. 

In Fig. HI we compare the intramolecular contact maps of A/540 and A/542 in the presence 
and absence of Els. Fig. H] shows the region containing the strongest contact V36-V39 as 
reported in our previous study (rectangle 1 in (a) and (c)) 32|. In A/540 the contacts 
between the amino acid regions L34-V36 and V39-V40 are significantly weaker than similar 
contacts between L34-V36 and V39-A42 in A/542. Els do not affect contacts in the TRB 
region (rectangle 1 in (b) and (d)). This result suggests that Els do not alter the contacts 
that contribute to differences between A/540 and A/542 folding in the CTR. 

A few important contacts in both alloforms in the MHR, concentrated around the 
strongest contact I31-L34, bring into proximity the two MHR regions A30-I32 and L34- 
V36 and are not affected by Els (rectangle 3 in (a)-(d)). The formation of these contacts 
within the MHR is promoted by G33 because glycines are typically associated with a high 
turn/loop propensity. Contacts between the CTR and MHR are present in both A/540 (rect- 
angle 2 in (a)) and A/542 (rectangle 2 in (c)), but are significantly stronger in A/542. These 
contacts are not affected significantly by Els (rectangle 2 in (b) and (d)). 

The central and most abundant contacts in folded monomers of both alloforms are formed 
as a consequence of the formation of the turn involving the TRA region (rectangles 4 and 7 
in (a)-(d)). The TRA region contains charged amino acids E22, D23, and K28, thus it is 
expected that Els will influence the contacts in this region. A strong contact A21-V24 in the 
TRA region becomes weaker as a result of Els (rectangle 7 in (a)-(d)), which is consistent 
with the effect of Els on the turn propensity in this region, changing a two-turn region into 
a one-turn region. Formation of contacts within the TRA brings into proximity the CHC 
and MHR (rectangles 5 in (a)-(d)). In both alloforms in the absence of Els, the CHC region 
makes contacts with the MHR with F19-I31 as the strongest contact (rectangles 5 in (a) 
and (c)). Els enhance the contacts within and around the TRA region in both alloforms, 
making contacts between the regions L17-D23 and K28-I32 (rectangles 5 in (b) and (d)) 
stronger. This enhanced feature is a consequence of a salt bridge formation between the op- 
positely charged D23 and K28. The TRA region was recently hypothesized to represent the 
nucleation re gion of A/5 folding 41] . This turn has been shown to be important in the fibril 



structure 



45 



461 ]. suggesting that this region maintains conformational stability throughout 



the folding and assembly of A/5. Our results are consistent with this hypothesis as they show 
that formation of contacts within the TRA region induces prominent contacts between the 
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CHC and MHR, resulting in the highest concentration of intramolecular contacts, involving 
the TRA, CHC, and MHR. 

In the absence of Els, the MHR region A30-M35 makes contacts with both the CHC 
(rectangle 5 in (a)-(d)) and CTR (rectangle 2 in (a)-(d)). These contacts do not change 
significantly in the presence of Els. The difference between A/340 and A/342 is that in A/340 
contacts between the regions A30-I32 and L34-V36 are stronger than the contacts between 
A30-I35 and V39-V40, while in A/342 the contacts between the regions A30-I35 and V39-A42 
are dominant. This result suggests that in A/342 folding the CTR plays a prominent role, 
while in A/340 the contacts within the MHR and between MHR and CHC regions are more 
important. 

The contacts between the K16-F19 and E11-H14 become more pronounced in the pres- 
ence of Els due to the EI between the negatively charged Ell and positively charged K16 
(rectangle 8 in (a)-(d)). A weaker group of contacts within the NTR between F4-H6 and 
Y10-V12 is a result of a turn centered at D7-G9 and hydrophobic attraction F4-V12. These 
contacts are very weak in the absence of Els (rectangle 9 in (a) and (c)) but become stronger 
in the presence of Els due to salt bridge R5-E11 (rectangle 9 in (b) and (d)). 

Long-range contacts between V39-V40 and CHC are observed in both A/340 and A/342 in 
the absence of Els (rectangle 6 in (a) and (c)). These contacts remain strong in the presence 
of Els (rectangle 6 in (b) and (d)). In A/342, these contacts are stronger than in A/340, 
both in the absence and presence of Els. Another region of long-range contacts is observed 
in both alloforms between the K16-F20 and D1-F4 in the absence of Els (rectangle 10 
in (a) and (c)). These contacts become more pronounced in the presence of Els due to 
electrostatic attraction between the negatively charged Dl and E3 and positively charged 
K16 (rectangle 10 in (b) and (d)). The long-range contacts between CTR and A2-F4, and 
MHR and A2-F4 are also present in both A/340 and A/342 but are weaker than the others 
and are not significantly influenced by Els. 

Time progression of A/3 folding events 

Fig. [5] shows time evolution of A/340 and A/342 monomer folding events in the presence 
of Els. Initially, A/340 and A/342 monomers are in zero potential energy, random coil con- 
formations. At lk simulation steps, contacts are formed between L34-V36 and CTR in both 
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A/340 and A/342. However, only in A/342, these contacts are associated with a turn structure 
in the TRB region as described in the previous section. At 2k steps, the contacts between 
regions CHC and TRA, CHC and MHR, CHC and CTR develop in both A/340 and A/342. 
These contacts are associated with a turn structure in the TRA region in both A/340 and 
A/942. At 4k steps, contacts between NTR and CHC develop in A/340. At 8k steps, as the 
contacts between NTR and CHC in A/340 are more pronounced, these contacts also emerge 
in A/342. At 0.1M steps, the long-range contacts between NTR and CTR are formed in both 
A/940 and A/342. Using the regions defined in Figs. |4] (b) and (d), the time progression of 
contacts follows the numbering 1, 2, 3, ... 10, i.e., A/3 folding starts at the C-terminal and 
progresses towards the N-terminal. In A/340, the turn structure in the TRA region is the 



first structura 
of Lazo et al. 



ement that is formed in the process of folding, supporting the hypothesis 
411 ] stating that the region 21-30 nucleates A/3folding. However, in A/342 the 
turn structure in the TRB region is formed before the formation of the turn structure in the 
TRA region. This result suggests that in A/342 the TRB region nucleates the folding prior 
to formation of contacts in the TRA region. 



Secondary structure of A/3 pentamers and larger oligomers 



In our previous work |32[], we reported the secondary structure difference between A/340 
and A/342 pentamers that can be found in the NTR and CTR. A/342 pentamers displayed an 
increased /3-strand propensity at the V39-I41, while A/340 pentamers showed an increased 
/3-strand propensity at the A2-F4. Our present data show that these differences remain 
intact in the presence of Els. 

Pentamers and larger oligomers do not have any significant amount of a-helix structure 
(data not shown). Figs. [6] (a) and (b) show the turn propensity per amino acid for A/340 
and A/342 pentamers and larger oligomers in the absence and presence of Els. Els do not 
affect the turn propensity significantly. In A/342, a slight increase in the turn propensity 
due to Els is found in the region R5-Y10. 

Figs. |6] (c) and (d) show the /3-strand propensity per amino acid for A/340 and A/342 
pentamers and larger oligomers. In both alloforms, the /3-strand propensity in the region 
K28-I31 slightly increases and in the region L34-G38 decreases due to Els. In the presence 
of Els, the /3-strand propensity in the CHC increases in A/340, while it decreases in A/342 
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pentamers and larger oligomers. 

We also calculate the average turn and /3-strand contents within A/540 and A/542 pen- 
tamers and larger oligomers in the absence and presence of Els. The data is shown in Table 
II. The average contents are calculated from propensities per residue by averaging over all 
residues in the peptide. The average turn content is in the range 41-45% and the average 
/3-strand content is in the range 11-13%. There is no significant difference between the two 
alloforms and no significant effect due to Els. 

These results show that pentamers and larger oligomers in our study have a globular 
structure dominated by turns and loop and some /3-strand propensity. Els change the 
relative /3-strand propensities of some regions, but do not affect significantly the overall 
secondary structure. 

Ramachandran plots of selected amino acids within the A/342 pentamers and higher 
oligomers 

Because our protein model as well as the interactions are simplified, we tested A/342 
oligomer conformers by calculating the Ramachandran plots. We selected the following 10 
amino acids from different regions of the protein: Dl, Y10, F19, E22, D23, S26, K28, M35, 
141, and A42. 

Our results shown in Fig. [7| indicate that both in the absence and presence of Els, the 
most populated ($, region corresponds to the /3-sheet region. The exceptions are Dl 
and A42, the N- and C-terminal amino acids, due to an increased flexibility at the two 
termini, and E22. Interestingly, E22 shows a substantially higher propensity to form a 
right-handed alpha-helix. Our results show that Els do not affect these plots in a significant 
way. These results are in qualitative agreement with A/3 dimer analysis of Huet et al. who 
studied A/3 dimer conformations by all-atom MD fl4j |. suggesting that our four-bead model 
yields relatively realistic set of <3> and \l/ angles and thus adequately accounts for the protein 
backbone structure. 
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Tertiary structure of pentamers and larger oligomers 



The tertiary structure of A/3 molecules within pentamers and larger oligomers (Fig. [8]) 
is highly reminiscent of the structure of individual monomers (compare Figs. H] and [SJ, sug- 
gesting that no major refolding events are needed in monomers prior to oligomer formation. 
However, there is less involvement of the N-terminal amino acids and more intramolecular 
contacts involving the C-terminal amino acids in A/3 molecules comprising pentamers and 
larger oligomers of both alloforms. 

There are significant differences between A/340 and A/342 intramolecular contact maps 
of pentamers and larger oligomers. The differences between A/940 and A/342 in the absence 



of Els have been described in our previous work [32| and can be observed comparing the 
relative importance of the CHC and CTR: in A/342 the contacts of CTR with MHR and 
CHC are dominant, while in A/340 the CHC plays a dominant role. In A/340 (Figs. E£a) 
and (b)) the contacts in regions marked by rectangles 1, 3, 4, and 5 get weaker due to 
Els, while the opposite is true in A/342 (Figs. [8](c) and (d)), where the contacts within the 
rectangles 1, 2, 3, 4, and 5 get stronger. This effect of Els on the intramolecular contacts 
can only be observed in pentamers and larger oligomers and not in unassembled monomers. 
A/342 pentamers and larger oligomers, in the presence of Els, have significantly stronger 
intramolecular contacts than A/340, suggesting that A/342 pentamers and larger oligomers 
are intrinsically more stable than their A/340 counterparts. 

Fig. [7] shows Ramachandran scattering plot on pentamers and larger oligomers of A/342. 
As seen from contact map analysis, in the presence of Els, Dls are more populated in /3-sheet 
region, which is the upper left corner. 



Quaternary structure of pentamers and larger oligomers 

Intermolecular contact maps indicate contacts among different A/3 molecules within an 
oligomer that are most important in oligomer formation. Previously, we showed that in 
A/340 pentamers, pairs of the CHC regions show the highest propensity to interact, whereas 
in A/342 pentamers the most frequent contacts are between the CTR of one peptide and 



the CHC and MHR of the other 



32( | . That result indicated that the CTRs are critically 



involved in aggregation of A/342 but not A/340. 
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Fig. [9] shows intermolecular contact maps of pentamers and larger oligomers of A/340 and 
A/542 in the absence ((a) and (c)) and presence ((b) and (d)) of Els. Perhaps the most 
surprising overall observation is that the intermolecular contacts that involve the CHC, i.e., 
contacts between pairs of CHCs (rectangle 3 in (a)-(d)), between the CHC and MHR (rect- 
angle 5 in (a)-(d)), and between the CHC and CTR (rectangle 6 in (a)-(d)), become weaker 
as a consequence of Els in both alloforms, but this weakening is more pronounced in A/340 
oligomers. This weakening of the contacts involving the CHC due to Els is surprising because 
the CHC is surrounded by charged residues (K16, E22, and D23). Thus, we would expect 
CHCs to interact pairwise in an anti-parallel fashion to maximize the the mutual attraction 
involving hydrophobic residues by additional salt bridge formation and thus minimize the 
free energy. Instead, our results show that Els weaken the contacts between pairs of CHCs. 
We also showed that Els promote formation of larger oligomers in both A/340 and A/542. 
These two results combined imply that weakening of the contacts between pairs of CHCs in 
A/340 oligomers might actually indirectly promote aggregation into larger oligomers. 

The only exception to the above observation is the region between D1-R5 and K16-D23 
which is rather weak in both alloforms in the absence of Els, but gets more pronounced in 
particular in A/342 due to Els (rectangle 7 in (a)-(d)). 

Our results indicate important differences in the way Els affect A/340 and A/342 oligomers. 
In A/340 oligomers the intermolecular contacts between pairs of CTRs (rectangle 1 in (a) 
and (b)), between pairs of MHRs (rectangle 2 in (a) and (b)), and between the CTR and 
MHR (rectangle 4 in (a) and (b)) remain unaffected by Els. In A/342 oligomers, on the other 
hand, the intermolecular contacts in these same regions get stronger even though that part 
of A/342 (MHR and CTR) is free of charge and thus Els would not be expected to make a 
difference. The strongest increase in the intermolecular contact intensity in A/342 oligomers 
is between pairs of CTRs (rectangle 1 in (b) and (d)) and the second strongest is between the 
CTR and MHR (rectangle 4 in (b) and (d)). Thus, in A/342 oligomers the contacts involving 
the CHCs get weaker and the contacts involving the CTRs get stronger due to Els, resulting 
in a significantly larger oligomers. These results suggest that in A/342 the CTRs are most 
important for intermolecular assembly into pentamers and larger oligomers. The lack of 
strong intermolecular contacts involving CTRs in A/340 suggests that the CTRs are also 
the main source of the differences between A/340 and A/342 oligomer formation. Recently, 
the importance of the intermolecular CHC contacts in A/340 versus the intermolecular CTR 



17 



contacts in A/342 was observed experimentally by Maji et al. 62J, in agreement with our 
present in silico results, suggesting the biological relevance of our DMD approach which is 
able to capture the essential differences between A/340 and A/342 oligomerization. 



Intra and intermolecular hydrogen bonds in pentamers and larger oligomers 

Here we address the question of how much hydrogen bonds contribute to intra- and in- 
termolecular contacts in pentamers and larger oligomers. We first calculate the probabilities 
for forming an intra- or intermolecular hydrogen bond per amino acid. The amino acids that 
are most hydrogen bond active are shown in Tables III and IV. Our results show that even 
for the amino acids that are most likely to form hydrogen bonds, probabilities are smaller 
than 0.20. The sum of intra- and intermolecular probabilities per amino acid does not exceed 
0.30/0.40, which is consistent with the /3-strand propensity per amino acid (Fig. [6]). 

Fig. [TU1 shows the intramolecular hydrogen bond contacts of A/340 ((a) and (b)) and A/342 
((c) and (d)) pentamers and larger oligomers in the absence ((a) and (c)) and presence ((b) 
and (d)) of Els. These intramolecular hydrogen bond maps are normalized to the highest 
value of intramolecular hydrogen bond formation probability, which is < 0.09. The regions 
with the highest amount of hydrogen bonds can be found between the regions K16-V24 and 
K28-V40. In A/342 oligomers some additional hydrogen bonds are formed between the MHR 
and CTR and between the CHC and CTR. Els increase the hydrogen bond probabilities 
within the TRA region and between the CHC and MHR due to salt bridge D23-K28. This 
effect is more pronounced in A/340. Interestingly, the strongest intramolecular hydrogen 
bond occurs in A/342 oligomers between F4 and V12, possibly stabilized by proximity of 
oppositely charged R5 and Ell. Why this same hydrogen bond is missing in A/340 oligomers 
may be understood by observation that in A/340 the region A2-F4 forms a /3-strand that is 
in contact with the CHC and thus the charged NTR residues (E3 and R5) are interacting 
with the charged residues K16 and E22, preventing R5-E11 from interacting and breaking 
the F4-V12 hydrogen bond. 

The intermolecular hydrogen bonds of A/340 ((a) and (b)) and A/342 ((c) and (d)) pen- 
tamers and larger oligomers in the absence ((a) and (c)) and presence ((b) and (d)) of Els 
are presented in Fig. [TTJ These intermolecular contact maps are normalized to the highest 
value of intermolecular hydrogen bond probability, which is < 0.04. The probability of in- 
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termolecular hydrogen bond formation is slightly higher in the regions where the contacts 
are more pronounced. Els do not influence the intermolecular hydrogen bond formation in 
any significant way. 

Our results show that the hydrogen bonds present in A/3 pentamers and larger oligomers 
are not specific, indicating that oligomers are not characterized by any particular pattern 
of hydrogen bonding. These findings suggest that hydrogen bonding is mostly a secondary 
effect occurring as a consequence of hydrophobic contact formation in the regions CHC, 
MHR, and CTR. 



Conclusions 



Because molecular dynamics approach to study proteins using all-atom representation 
and explicit solvent is limited to time scales smaller than ~ 10~ 6 s, we use a simplified but 
efficient DMD approach combined with a four-bead protein model and amino acid-specific 
interactions that mimic the effects of a solvent llj . In our prior work we showed that this 
approach yields biologically relevant results, which are consistent with existing experimental 
findings on A/3 oligomer formation and have predictive power allowing for in vitro and further 



in silico testing [32|. In the present work we use the DMD approach to study the effects 
of Els on oligomer formation of A/540 and A/342. The role of electrostatic interactions, 
in particular the salt bridge formation between negatively charged E22/D23 and positively 
charged K28 was hypothesized to be important at early stages of folding as well as at later 
stages of fibril formation. Thus, it is reasonable to expect that Els may play an important 
role at intermediate stages of oligomer formation. 

We analyze the structure of folded A/340 and A/342 monomers in the presence and absence 
of Els. We show that independent of Els the two alloforms display differences in folded 
structure: in A/342 there is an additional turn centered at G37-G38 that is absent in A/340, 
leading to an increased propensity to form /3-strand in the CTR of only A/342. A/340 
monomers also have an additional /3-strand in the A2-F4 which is not present in A/342. Our 
results demonstrate that the differences between the two alloforms are present already at 
the stage of folding, prior to assembly. The existence of a turn structure centered at G37- 
G38 is consistent with experimental findings by Lazo et al. who showed by using limited 
proteolysis that Val39-Val40 in A/342 but not in A/340 monomer was protease resistant, 
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indicating that A/342 but not A/340 monomer was structured in the CTR region 41J . Similar 
was a conclusion of the solution NMR study on [Met(O) 35 ]A/340 versus [Met(0) 35 ]A/342 
monomer structure by Riek et al. showing that G29-A42 region is less flexible and thus 



more structured in A/342 than in A/340 



631 ] . By measuring 1 H a , l3 C a , and Vi Cp chemical 



64j which provides further 



shift indices of A/340 and A/342 Hou et al. recently showed that the C-terminus of A/342 but 
not of A/340 monomer has a tendency to form /3-sheet structure 
evidence that our simulation approach yields biologically relevant results consistent with in 
vitro findings. 

Our results indicate that Els stabilize a turn in the region D23-K28 by formation of a 



D23-K28 salt bridge. A role for Els in stabilizing this region has been postulated 



et al. 



jy Lazo 



411 ] and further explored using a more complex united atom DMD model 



all-atom MD in explicit 



43j and 



44j and implicit solvent 29(. These studies show that A/3 folding 



in the region A21-A30 is driven primarily by effective hydrophobic attraction between V24 
and the butyl portion of K28, but that Els help stabilize the region. In our model, due to its 
simplicity, the side chains of V24 and K28 do not experience attractive interactions. Despite 
the absence of this interaction, we still find this region to be the most structured in both 
A/340 and A/342 monomers stabilized by D23-K28 salt bridge. The D23-K28 salt bridge was 
suggested to stabilize the A/340 fibril structure by Petkova et al. 46j . In addition, Sciarretta 
et al. have shown an increase in the rate of A/MO-Lactam (D23/K28) fibrillogenesis by 1000- 
folds (47)], providing additional experimental evidence supporting a critical role of D23-K28 
salt bridge formation. 

Comparing the oligomer size distributions of A/340 and A/342 in the presence of Els 
with those obtained in the absence of Els 32J reveals that Els promote formation of larger 



oligomers while maintaining a unimodal A/340 size distribution and a multimodal A/342 
size distribution, as observed in vitro In our simulations the N- and C-termini are 
uncharged in contrast to most experimental studies with positively charged N- and negatively 
charged C-termini. Our observation that Els promote formation of larger oligomers is thus 
consistent with results of the experimental study in which the C-terminal carboxyl group 
was replaced by the electrostatically neutral carboxamide, resulting in a greater abundance 
of high molecular weight oligomers 42J . 

It is critical to study the structural changes in oligomers due to Els and understand which 
structural changes are contributing to formation of larger oligomers in both A/340 and A/342. 
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Our results indicate that in A/340 pentamers and larger oligomers, Els weaken intramolecular 
interactions. In A/342, in contrast, the intramolecular contacts in the turn region D23-K28 
are enhanced. Surprisingly, in both A/340 and A/342 oligomers, the intermolecular contacts 
involving the CHC are significantly weaker in the presence of Els. In addition, in A/342 
oligomers, the contacts involving the CTR and MHR get stronger. These results, combined 
with the fact that Els promote larger oligomers, imply that the intermolecular interactions 
between pairs of CHCs in an indirect way oppose the formation of larger oligomers, while 
the interactions between pairs of CTRs, and to a smaller extent also pairs of MHRs, promote 
formation of larger oligomers. Thus, therapeutic strategies using inhibitors that target the 
CTR and MHR may prove successful in either inhibiting formation of toxic A/342 oligomers 
or inducing structural modifications neutralizing their toxicity. 
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Figure captions: 

Figure 1: Representative conformations of a monomer, dimer, and pentamer of A/542 in 
the absence (Ech — 0) and presence (Ech = 0.6) of Els. A monomer conformation (a) 
in the absence and (b) presence of Els. A dimer conformation (c) in the absence and (d) 
presence of Els. A pentamer conformation (e) in the absence and (f) presence of Els. Yellow 
arrows correspond to the /3-strand structure, turns are represented by light blue tubes and 
random coil-like part are represented by gray tubes. The N-terminal Dl is marked as a 
red sphere, and the C-terminal A42 is marked as a blue sphere. 131, 132, and 141, the most 
hydrophobic residues, are represented as green spheres. This figure is generated by the VMD 
package [65|. 

Figure 2: Oligomer size distributions of A/340 and A/342 at (a) Ecu — and (b) Ech — 
0.6. All size distributions are averages over the time frames at 8, 8.5, 9, 9.5, and 10 million 
simulation steps. 

Figure 3: The effect of Els on turn and /3-strand propensities per residue in folded A/3 
monomers in the absence and presence of Els. Turn propensities of (a) A/540 and (b) A/342 
monomers. /3-strand propensities of (c) A/340 and (d) A/342 monomers. 

Figure 4: Intramolecular contact maps of folded A/340 and A/342 monomers at Ecu — 
(left column) and Ech = 0.6 (right column). The strength of the contact map is color- 
coded following the rainbow scheme: from blue (no contacts) to red (the largest number 
of contacts). Each contact map is an average of over 100 monomer conformations after 1 
million simulation steps. 

Figure 5: Detailed time evolution of intramolecular contact formation during A/340 (left 
column) and A/342 folding (right column). The strength of the contact map is color-coded 
as in Fig. 4. Each contact map is an average of over 100 monomer conformations. 

Figure 6: The effect of Els on turn and /3-strand propensities per residue within A/3 
pentamers and larger oligomers in the absence and presence of Els. Turn propensities of (a) 
A/340 and (b) A/342 pentamers and larger oligomers. /3-strand propensities of (c) A/340 and 
(d) A/342 pentamers and larger oligomers. 

Figure 7: Ramachandran plots of A/342 pentamers and larger oligomers for selected 
residues Dl, Y10, F19, E22, D23, S26, K28, M35, 141, and A42 in the absence and presence 
of Els. Horizontal and vertical axes correspond to the angles $ and respectively, both 
varying from -180° to 180°. Each plot contains ~640 points corresponding to A/342 pen- 
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lion simulation steps. Ramachandran 
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tamers to decamers obtained at 8, 8.5, 9, 9.5, and 10 mi 
plots are generated using the VMD software package 

Figure 8: Intramolecular contact maps of A/340 and A/342 pentamers and larger oligomers 
at Ecu — and Ech — 0.6. Each contact map is an average of over 50 conformations 
obtained at 8, 8.5, 9, 9.5, and 10 million simulation steps. 

Figure 9: Intermolecular contact maps of A/340 and A/342 pentamers and larger oligomers 
at Ech = and Ech = 0.6. Each contact map is an average of over 50 conformations at 8, 
8.5, 9, 9.5, and 10 million simulation steps. 

Figure 10: Intramolecular hydrogen bond maps of A/340 and A/342 pentamers and larger 
oligomers at Ech = and Eqh — 0.6. Each map is an average of over 50 conformations at 
8, 8.5, 9, 9.5, and 10 million simulation steps. 

Figure 11: Intermolecular hydrogen bond maps of A/340 and A/342 pentamers and larger 
oligomers at Ech — and Ech — 0.6. Each map is an average of over 50 conformations at 
8, 8.5, 9, 9.5, and 10 million simulation steps. 

Table captions: 

Table 1: Average turn and /3-strand propensities per residue with standard errors within 
folded A/340 and A/342 monomers. Each value is an average of over 100 monomer confor- 
mations after 1 million simulation steps. 

Table 2: Average turn and /3-strand propensities per residue with standard errors within 
A/340 and A/342 pentamers and larger oligomers. Each value is an average of over 50 con- 
formations at 8, 8.5, 9, 9.5, and 10 million simulation steps. 

Table 3: Average hydrogen bond propensities per residue, showing the five most frequent 
residues involved in intramolecular hydrogen bonding within A/340 and A/342 pentamers and 
larger oligomers. Each value is an average of over 50 conformations at 8, 8.5, 9, 9.5, and 10 
million simulation steps. 

Table 4: Average hydrogen bond propensity per residue, showing the five most frequent 
amino acids involved in intermolecular hydrogen bonding within A/340 and A/342 pentamers 
and larger oligomers. Each value is an average of over 50 conformations at 8, 8.5, 9, 9.5, 
and 10 million simulation steps. 
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40 





A/540 


A/342 




Ech = 


Ech = 0.6 


E CH = 


Ech = 0.6 


Turn 


0.44 ± 0.04 


0.43 ± 0.04 


0.48 ± 0.04 


0.50 ± 0.05 


/?-strand 


0.11 ± 0.02 


0.12 ± 0.03 


0.11 ± 0.03 


0.10 ± 0.03 
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A/540 


A/342 




Ech = 


Ech = 0.6 


E CH = 


Ech = 0.6 


Turn 


0.44 ± 0.02 


0.45 ± 0.02 


0.42 ± 0.02 


0.45 ± 0.02 


/?-strand 


0.13 ± 0.02 


0.11 ± 0.01 


0.13 ± 0.01 


0.11 ± 0.01 
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A/340 


A/542 


E C h = 


Ecu = 0.6 


Ecu = 


Ecu = 0.6 


L17 


0.17 


M35 


0.14 


V40 


0.17 


A30 


0.17 


A21 


0.15 


G38 


0.14 


131 


0.16 


G29 


0.15 


G33 


0.15 


131 


0.13 


G38 


0.15 


Ell 


0.15 


V24 


0.14 


G33 


0.12 


A21 


0.13 


R5 


0.14 


G38 


0.14 


G37 


0.12 


G29 


0.13 


131 


0.14 
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A/340 


A/542 


E CH = 


Ecu = 0.6 


E C h = 


E C h = 0.6 


M35 


0.15 


M35 


0.14 


F20 


0.18 


131 


0.14 


131 


0.14 


G38 


0.14 


V18 


0.15 


G33 


0.13 


V36 


0.14 


131 


0.13 


G33 


0.13 


V40 


0.13 


G37 


0.14 


G33 


0.122 


141 


0.13 


L17 


0.12 


G38 


0.13 


G37 


0.12 


A21 


0.12 


F20 


0.12 
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